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Y. Dubi and M. Di Ventra 
Department of Physics, University of California San Diego, La Jolla, California 92093-0319, 



USA 



Despite its intrinsic non-equilibrium origin, thermoelectricity in nanoscale systems is usually de- 
scribed within a static scattering approach which disregards the dynamical interaction with the 
thermal baths that maintain energy flow. Using the theory of open quantum systems we show 
instead that unexpected properties, such as a resonant structure and large sign sensitivity, emerge 
if the non-equilibrium nature of this problem is considered. Our approach also allows us to define 
and study a local temperature, which shows hot spots and oscillations along the system according 
to the coupling of the latter to the electrodes. This demonstrates that Fourier's law - a paradigm 
of statistical mechanics - is generally violated in nanoscale junctions. 
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Non-equilibrium (NE) processes at the nanoscale are 
receiving a great deal of attention due in large part to the 
advancements in fabrication and manipulation of these 
systems. [1] An especially interesting class of NE phe- 
nomena pertain to energy transport and the conversion 
of thermal to electrical energy. When a thermal gradient 
AT is applied to a finite system, electrons respond by 
departing from their ground state to partially accumu- 
late at one end of the system, thus creating a measurable 
voltage difference AV . The ratio S = — ^ is called 
thermopower [2], and has been measured in a variety 
of nano-scale systems such as quantum point contacts 
[3], atomic-size metallic wires [4], quantum dots [5], Si 
nanowires [6] and recently in molecular junctions [7]. In 
a bulk material, when S < the transient current is 
carried by electrons; when S > it is carried by holes. 

In nanoscale systems this NE problem has recently re- 
ceived a lot of attention [3, 8-13]. In these theories the 
single-particle scattering formalism [14] is used to relate 
the thermopower to single-particle transmission proba- 
bilities. This approach, however, does not take into ac- 
count the dynamical formation of the thermopower and 
neglects the fact that even at steady state, when the 
charge current is zero an energy current is still present, 
like, e.g., in insulators [15]. Another effect neglected by 
such theories, which is now within reach of experimen- 
tal verification [16], is the formation of local temperature 
variations along the structure. In order to study all these 
effects one needs to describe a nanoscale system inter- 
acting with an environment that maintains the thermal 
gradient, namely one needs to resort to a theory of NE 
open quantum systems. 

In this letter we introduce such a theory, based on a 
generalization of quantum master equations, and use it to 
study the dynamical formation of thermo-electric effects 
in nano junctions. We show that the thermopower is a 
highly non-linear function of the thermal gradient and 
it is very sensitive to the junction geometry, even in the 
simplest case of non- interacting electrons. This precludes 
an easy interpretation of its sign in terms of electrons or 
holes as it has been argued in some literature [3, 8-11]. 



In addition, we calculate the global and local electron 
distribution functions, which exhibit NE characteristics. 

The theory also allows us to define the local electron 
temperature by means of a temperature floating probe 
that is locally coupled to the system, and whose tempera- 
ture is adjusted so that the system dynamics is minimally 
perturbed. This temperature, which can be measured ex- 
perimentally, shows important features such as hot spots 
in the cold lead at small coupling between the nanowire 
and the bulk electrodes, and temperature oscillations in 
the wire at intermediate coupling. These findings show 
that Fourier's law, which is considered a paradigm of 
thermodynamics, is generally violated for electronic sys- 
tems at the nanoscale [17]. 

Method. - Since we consider non-interacting electrons 
coupled to an environment, we employ a quantum mas- 
ter equation of the Lindblad type [18] which describes 
the dynamical evolution of the many-body density ma- 
trix (DM) pM of a quantum system in the presence of a 
markovian bath, via the introduction of a super-operator 
C[pm] [19]- The quantum master equation is then (h = 1) 



Pm —i[H, Pm] + C[p M ] 



(1) 



where [•, •] denotes the commutator. The super-operator 
C is defined via a set V nn > of operators via 

C\PM] = (-^{VL'Vnn>,PM} + V nn , PM V^ n ) , (2) 

n,n' 

with {•, •} being the anti-commutator. The sums over 
n and n' (n ^ n') are performed over all many-particle 
levels of the system, and the F-operators are conveniently 
selected in the form V nn > = ^^ nn ' \^ n )(®n' I? describing a 
transition from the many-body state \^ n >) into the state 
\Sfr n ) with the transition rate 7 rm / - 

This problem scales exponentially with the number 
of particles, but recently a mapping of the many-body 
super-operator to a single-particle form has been intro- 
duced [20]. This results in a quantum master equation 
for the single- particle DM, p = J2 kk , pkk>\k>)(k'\ which 
provides excellent agreement with the many-body solu- 
tion. Here, \k) are the single-particle states, and the 
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matrix elements are derived from the many-body DM by 
Pkk' = Tr (c\c kl pM^j. 

To be specific, we consider a finite nano-junction (i.e., 
with a fixed number of electrons and ions) which is com- 
posed of two identical quasi-two-dimensional leads con- 
nected via a one-dimensional wire (see upper panel of 
Fig. 1). The far edges of the leads are coupled to two dif- 
ferent baths kept at different temperatures. The Hamil- 
tonian of the system is given by TL = 'HL + TLR + TLd + 'Hci 

where TLl,R4 = ~^<;,j)eL,i?,d ( c l c i + h - c -) are tne 
tight-binding Hamiltonians of the left lead, right lead and 
wire, respectively (t is the hopping integral, which serves 
as the energy scale hereafter), and H c = (gLC L Cd,o + 
9R,c\iCd,L d + h.c.) describes the coupling between the left 
(right) lead to the wire, with being the creation 

operator for an electron at the point of contact between 
the left (right) lead and the wire, and q,o ( c d,L d ) destroys 
an electron at the left-most (right-most) sites of the wire. 
We consider here spinless electrons. The master equation 
now takes the form 



p=-i[H,p] + £ L \p] + CR\p] 



(3) 



where jCl(r) describes relaxation processes due to the 
contact between the left (right) lead with its respective 
bath at temperature T L ^ R y The V-operators are gener- 
alized to account for the different baths, and are given 
by [20, 21] 



(l,r)Al,r), 



(4) 



where f£" R) (e k ) = 1/ (exp ( fc ^VA ) + X ) are the Fermi 
distributions of the left and right leads, with /i the chem- 
ical potential. The coefficients 



(L,R) 
^kk' 



^fc( r i)7o^jfe/(ri) 



(5) 



describe the overlap between the single-particle states \k) 
and \k') over the region of contact S^i?) between the 
left (right) baths and the corresponding junction leads, 
shown by the solid lines in the upper panel of Fig. 1. 70 
describes the strength of electron-phonon (bath) interac- 
tion. The form (5) can be derived from first principles by 
tracing out the bath degrees of freedom, with the latter 
formed by a dense spectrum of boson excitations (e.g., 
phonons), which interact locally with electrons at the 
edges of the system. The operators (4) guarantee that 
the system evolves to a global equilibrium if Tl = Tr, 
or equilibrate each lead at its own temperature if g = 0, 
i.e., the leads and the wire are completely decoupled (and 
hence no voltage drop can form). 

We now solve equation (3) numerically for several tem- 
perature gradients. From the obtained charge density 



distribution we derive the electrical potential via the 
Poisson equation. The potential is averaged along the 
transverse direction, and the voltage drop is calculated 
from the center of the leads [22]. The off-diagonal ele- 
ments of the DM decay fast (on a time-scale ~ j^ 1 ~ 10 
in our calculations) and hence do not contribute to the 
density in the long-time limit. This allows us to neglect 
them completely, a fact which significantly simplifies the 
calculation [23]. 

Numerical results. - In Fig. 1(a) the voltage drop AV 
across the junction is plotted as a function of the temper- 
ature difference AT between the contacts. The leads are 
of dimensions 12 x 11 and the wire is of length Ld = 6. 
The lead- wire coupling is = qr = 0.001 and the num- 
ber of electrons is tie = 90, which corresponds to | fill- 
ing. The initial temperatures are set to Tl = Tr = 0.05. 
From Fig. 1 one notices three regimes in the range of AT. 
At small AT, a linear-response regime can be identified. 
This is followed by a regime of rapid rise in AV, eventu- 
ally reaching a saturation at large AT, due to the finite 
size of the system. The solid line is a fit to an expo- 
nential rise. Although the parameters of the exponential 
fit depend on sample parameters, we found that the ex- 
ponential form is an excellent fit for all non-interacting 
junctions. In the inset of Fig. 1(a) we plot the generalized 
thermopower, S = — ^( at) , which reduces to the regular 
thermopower in the linear regime. As seen, S exhibits 
a resonance at AT « 0.25 (this value is not universal 
and depends on junction parameters), which means that 
at this value the response of the system to a change in 
the temperature gradient is maximal; a fact that can be 
checked experimentally and may be used in actual de- 
vices. 

Geometrical effects. Due to the local variations of 
the density at the junction and hence local variations 
of kinetic energy, the thermo-electric response strongly 
depends on junction geometry, as it was anticipated ex- 
perimentally [4]. As an example, we have calculated the 
charge-imbalance AQ, across a junction (leads size 6x7, 
wire length Ld — 6, density at \ filling) with an asym- 
metric coupling between the leads and the wire. The 
coupling to the left lead was kept at gL = 0.001 and the 
coupling to the right lead, gn, was changed. In Fig. 1(b) 
we plot the charge imbalance as a function of gR at a 
fixed temperature difference AT = 1. Strong and nar- 
row oscillations can be seen, and for certain values of 
coupling asymmetry, AQ may even change sign. This is 
consistent with the experiments in [4] and may account 
for the sign change of the thermo-power observed in some 
experiments in molecular junctions [4, 7]. 

Nonequilibrium distributions. - The formalism presented 
here allows us to calculate various NE properties. As 
an example, we calculate the local and global distribu- 
tion function (DF) of the wire (which are accessible ex- 
perimentally [24]), for a system with lead size 5x5, 
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FIG. 1: (Color online) (a) Upper panel: nanojunction geom- 
etry considered in this calculation (see text for parameters). 
Main panel: Voltage drop AV across the junction as a func- 
tion of the temperature difference between the leads. There 
are three distinct regimes: the linear regime, a regime of rapid 
voltage rise and saturation. Solid line is a fit to an expo- 
nential rise. Inset: generalized thermopower, S = — ^[at) • 
(b) Charge imbalance AQ as a function of the ratio between 
the right and left lead- wire coupling gii/gL, at a constant 
QL 0.001 and at AT = 1. 



wire length = 160 (at third filling), with coupling 
9l = 9r = 1, and temperatures Tl = 0.05 and Tr = 0.5. 
In Fig. 2(a) we plot the full distribution function of 
the wire (circles), f(Ek) = pkk, along with the curve 
f(E k ) = \tf D {T L ) + f D {T R )) (solid line). As seen, the fit 
between the data and the DF is excellent. In Fig. 2(b-f) 
we plot the local DF f(Ek,x) = ptk 1^(^)1 2 f° r different 
positions along the wire, x = 1,40,80,120,160 (respec- 
tively). Note the fast oscillations upon approaching the 
center of the wire, and the symmetry with respect to the 
wire center. The origin of the oscillations lies both in the 
NE nature of the system and in the geometry (i.e. the 
fact that the wire is ballistic and that the wire-lead cou- 
pling is large). In a system with small wire- lead coupling 
we found (not shown) that the oscillations persist at the 
center of the wire, but are smoothed at the wire edges. 
Fig. 2(b-f) shows that although the full DF is a simple 
average of the left and right lead distributions, one can- 
not assign a simple position dependence, as in the case 
of mesoscopic wires [24] . 

Local temperature. The concept of a local temper- 




FIG. 2: (Color online) (a) Full distribution function of a long 
wire, Ld — 160 (see text for numerical parameters). The 
solid line is a mean distribution function f(E) — |(/d(Tl) + 
/d(Tr)). (b-f) Local distribution function for different posi- 
tions along the wire. 



ature is generally not unique out of equilibrium [25]. 
Here we provide an operational definition, which is both 
physically transparent and can (in principle) be directly 
probed experimentally [16]. In order to do so, we add 
an additional relaxation operator Aip [pm] to the master 
equation (1). This corresponds to applying a local bath 
at a temperature T t i p in contact with a single site of the 
system (see upper panel of Fig. 3). Due to energy flow 
between the probe and the coupled site, the system dy- 
namics is generally modified. However, one can "float" 
the temperature T t i p such that the change in local (and 
hence global, e.g., thermopower) properties of the system 
is minimal. We define this temperature as the local tem- 
perature of our system at the probe position. We choose 
to monitor the change in the local density from its value 
in the absence of the probe, but any other quantity would 
be equally valid and lead to the same local temperature. 

Knowledge of the local temperature allows us to ad- 
dress the problem of validity of Fourier's law at the nano- 
scale. It has been conjectured [17] and demonstrated for 
systems such as spin-chains [21, 26] and coupled har- 
monic oscillators [27] that in ballistic quantum systems 
Fourier's law is invalid. Stimulated by this open prob- 
lem, we plot in Fig. 3 the local temperature at steady 
state for three different values of the lead-wire coupling, 
g = 0.001,0.1,0.8 and a junction with lead dimensions 
4x3, wire length Ld = 21, temperatures Tl = 0.05 and 
Tr = 1.5, and electron density at third filling. For weak 
coupling (g = 0.001), we find that the temperature in- 
side the wire is very low, but a "hot spot" develops in 
the cold lead. As the coupling increases the hot spot van- 
ishes, and temperature oscillations develop in the wire. 
As expected, at high coupling, the wire equilibrates at 
a temperature that is roughly the average between Tl 
and Tr. For large lead- wire coupling the temperature in 
the wire is uniform, and most of the temperature drop 
occurs at the contacts, similar to what has been argued 
for the phenomenon of local ionic heating [28]. We note 
that when the coupling g ~ 1 , so that the wave-functions 
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FIG. 3: (Color online) Upper panel: schematic representa- 
tion of the calculation to determine the local temperature via 
the addition of a local temperature floating probe. Lower 
panels: the local temperature along the structure for three 
different values of the lead-wire coupling, g = 0.001,0.1,0.8. 
Three effects are observed: a hot-spot in the cold lead at small 
coupling, temperature oscillations in the wire at intermediate 
coupling, and uniform temperature along the wire at large 
coupling. All three constitute violations of Fourier's law. 



are completely delocalized along the nano-j unction, the 
thermopower is small. This is consistent with the exper- 
imental results of [4], where junctions of large conduc- 
tance show small thermopower values. Our results verify 
that Fourier's law is generally violated in nanoscale junc- 
tions. 

Summary. - We have shown, using an open quantum 
system approach, that the proper description of energy 
flow in nanoscale systems leads to unexpected features in 
the dynamical formation of thermopower. Several pre- 
dictions have been made which can be tested experimen- 
tally, namely (i) the non-linear dependence (with reso- 
nant structures) of the thermopower as a function of the 
temperature gradient, (ii) its strong sensitivity (includ- 
ing sign) to junction geometry, (iii) the shape of the NE 
electron distribution at steady state, and (iv) we have 
provided an operational definition of local temperature 
as that measured by a temperature floating probe, i.e., 
one that is locally coupled to the system, and whose tem- 
perature is adjusted so that the system dynamics is min- 
imally perturbed. This temperature shows noteworthy 
features according to the strength of the coupling be- 
tween the nanoscale system and the electrodes, in viola- 
tion of Fourier's law. 



We conclude by pointing that the results presented 
here may be relevant to other systems of present interest 
(e.g. graphene nano-ribbons, nanotubes etc.). Studies 
of these effects that include electron interactions repre- 
sent another important research direction, and are un- 
derway using stochastic time-dependent current density- 
functional theory [29]. 
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